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ABSTEIACT 
In  the  problem  of  evsLLuatlng  two  and.  three-center  integrals  which 
give  rise  to  infinite  series,  the  one  and  two-electron  integrals  (\rtiich 
are  three  and  six-dimensional  integrals)  are  reduced  to  infinite  series 
of  one  and  two-dimensional  integrals  respectively.  Convergence  rates 
euid  asymptotic  behavior  for  large  order  terms  of  the  series  are  derived. 
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1.   Introduction 

The  expression  of  the  solutions  of  SchrSdinger 's  equation  for  the 
electronic  configuration  of  a  molecule  as  a  linear  ccMnbination  of  given 
atomic  orbital  functions  transforms  Schrbdinger ' s  equation  into  a  matrix 
equation.  Tbe   solution  of  this  matrix  equation  involves  the  evaluation 
of  certain  three  and  six-dimensional  integrals  referred  to  as  one  and 
two- electron  IntegraLLs,  respectively.  Ihe  domains  of  integration  of 
the  integrals  are  all  of  Euclidean  three  emd  six  space.  Scane  of  the 
general  techniques  for  evaluating  these  Integrals  involve  integrating 
over  certain  of  the  variables,  transforming  the  integrals  into  series 
of  one  and  two-dimensional  integrals. 

The  exact  fonn  of  these  series  depends  upon  the  choice  of  coordinate 
system  in  Euclidean  three  space.  Once  the  choice  of  coordinate  system 
has  been  made  and  the  integrations  carried  out,  the  integrals  fall  into 
two  categories:  those  involving  infinite  series  and  those  involving 
finite  series. 

Two  methods  eire  camnonly  used  In  the  evaluation  of  these  integrals, 
each  based  upon  a  particular  choice  of  coordinate  system  in  Euclidean 
three  space.  One  method  involves  carrying  out  the  Integrations  using 
prolate  spheroidal  coordinates.  The  evaluation  of  the  nonocenter  and 
two-center  Integrals  using  this  method  has  been  carried  out  by  some 
earlier  workers  t^'^' 5'^' '''■1 . 

IHie  other  method  involves  carrying  out  the  integrations,  using 
spherical  coordinates.  The   evaluation  of  the  two  center  integrsils 
involving  a  finite  series,  (the  one  electron,  two  electron  monocenter, 
and  two  center  hybrid  and  coulosnb  integrals),  using  this  method  has  been 
carried  out  by  Bamett  and  Coulson  ^  -^ . 
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Bie  advantages  of  the  first  method  over  the  second  In  the  evaluation 
of  these  integrals  are:  (l)  more  of  the  two  center  integrals  involve  only 
a  finite  series,  and  (2)  the  Infinite  series  obtained  appear  to  converge 
rapidly.  Ihe  advantages  of  the  second  method  over  the  first  are:(l)  the 
special  functions  are  relatively  easy  to  coinpute,  and  (2)  the  second 
method  can  more  readily  be  extended  to  the  evaluation  of  three  and  four 
center  integrals  (arising  from  problems  in  the  molecular  structure  of 
molecules  with  three  and  four  atoms ) . 

In  this  paper  we  shall  be  concerned  with  the  evaluation  of  some  of 
the  two  and  three  center  integrals  of  molecular  structxire.  Because  we 
wish  to  evaluate  all  of  the  three  center  integrals,  we  shall  use  the 
method  involving  spherical  coordinates.  Using  this  method  we  shall 
evaluate  those  two  aiwi  three  center  integrals  which  give  rise  to 
infinite  series,  namely: 

(1.1)  J  t^"'"  (l)rci  *B  (1)  dv  (potential  energy  integrals) 

(1.2)  J  j  i^^)(l)Mr^2)(i)r-l  ♦p\2)*^'*^2)dv^dv2 
2  1 

(two  center  exchange  integrals) 

(1.3)  jj  tS,^\l)*^^^l)r^*P^(2)>lr^^)(2)dv^dv2 
2'1 

(three  center  hydrid  exchange  integrals) 

(l.lf)   JJ  ♦^^^l)*j^2)(^j^-1^0)(2)^(^l^)(2)dv^dv2 
2  1 

(three  center  coulomb  exchange  integrals) 


P3=  AB 


<  CAB  =  a 


PC=  AC 


FIGURE  1 


In  order  to  define  precisely  what  these  integrals  are,  ve  may 
consider  the  case  of  a  general  triatomic  molecule  in  three  dimensional 
space  (see  Figure  l),  where  A,  B,  C  are  three  arbitrary  noncolinear 
points  in  space  (the  locations  of  the  tliree  nuclei  of  the  atoms  of 
the  molecule) .  We  define  three  rectangular  coordinate  systems 
(x  ,  y  ,  z  )  ^  »  A,  B,  or  C  with  A,  B,  C  as  origins  as  follows. 
y«^  yjif  Jo   iiave  the  same  direction  ea*d  are  perpendicular  to  the 

A    O    w 

plane  determined  by  A,  B,  C.  z  is  taken  along  AB  directed  toward  B 
and  z_,  Zp  are  taken  along  AB,  AC  respectively,  directed  toward  A. 
X.,  Xg,  x_  are  chosen  so  that  (x^,  y^,  z^)  is  a  right-handed  coordinate 
system  and  (x,y,z)n=B,  C  are  left-handed  coordinate  systans.  In 

r*     r*     r* 

terms  of  these  rectangular  coordinate  systans  we  may  define  spherical 
coordinate  systems  (r  ,  e  ,  0  )  with  origins  at  ji  »  A,  B,  C  respectively. 
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with  these  definitions,  then,  we  may  define  the  integrals (1.3)  - 
(1.4).  The  +'^^(5);  1  «=  1,2,3,^}   n  «  A,B,C;  6  =  1,2  are  among  the  Slater 
functions  given  in  Appendix  (l) .  Here  we  restrict  ourselves  to  Slater 
functions  involving  major  quantum  numbers  n  =  1  or  2.   (^  5^  ©  =*  0  «) 
\i  "  A,   B,  C;  6  =  1,2  are  the  spherical  coordinates  of  electron  6  with 
origin  ii.     r.   is  the  distance  of  electron  2  from  electron  1  and  an 
expansion  for  r"  is  given  in  Appendix  (l).  dv  ,  5  =  1,2  sire  volume 
elements . 

2.   Method  of  evaluation  of  the  integrals 

The  method  of  evaliiating  the  integrals  used  in  this  paper  is  usually 
referred  to  as  the  Bessel  Function  Method  l-  -",  and  consista  essentially 
in  the  expansion  of  functions  with  argxnnents  (r  ,  9  ,  0  )  |.i  -  B  or  C 
in  terms  of  functions  with  arguments  ( r  ,  9  ,  0  ) .  The  expansions  used 
for  this  purpose  are  given  in  Appendix  (2),  and  their  derivation  Is 
straightforward • 

The  use  of  the  Bessel  Function  Method  in  evaluating  Integrals (1.1) - 
(l»  gives  rise  to  an  infinite  series  of  single  integrals  in  the  case  of 
integral  (l.l)  and  infinite  series  of  double  integrals  in  the  case  of 
integrals  (l.2)-(l.4).   It  Is  the  purpose  of  this  paper  to  establish 
the  rates  of  convergence  of  these  series  and  to  obtain  the  asymptotic 
behavior  of  the  general  terms  of  large  order  in  these  series. 

The  actual  techniques  for  reducing  the  three  euad  six  dimensional 
Integrals  to  an  infinite  series  of  one  and  two  dimensional  integrals 
respectively  and  the  analysis  of  the  asymptotic  behavior  and  the  con- 
vergence rates  of  the  series  eire  practically  identical  for  the  four 
cases,  but  are  somewhat  tedious  and  difficult  to  explain  without  using 
one  of  the  four  cases  as  an  example .  Therefore  the  general  techniques 
and  methods  will  be  explained.  Integral  (1.1)  as  an  example,  and  the  ex- 


pleination  will  not  be  repeated  in  the  other  three  cases  xuiless  a 
variation  in  the  techniques  Is  encountered. 

5.       Potential  energy  integrals 
"Hiese  are  of  the  form: 


(5.1)      |+^^^(l)r-^Mr^2)(i)dv^ 


=  C  C    X:       X:    E-    P^  (cos  a)  [    f    ( 


""ai      ® 

j  =  0  -I''  -1  ^Q  /q  ^Q 


^^/^f)  ^i,( 


Xr^fV'^'^^fp^/p.J     iPi(coBe^) 


X  P^     (cos  e^)  P^^     (cos  Gg^)  sin  e^ 

Starting  with  the  0.,  variable,  we  note  by  the  formulas  for  expanding 
products  of  sines  and  cosines  that  the  Integral  must  vanish  unless  the 
expansion  of  ^  (0  ,  )^  (0  )  io  terms  of  a  linear  coniblnatlon  of 
sines  eind  cosines  contains  a  term:  cos  M0^,  M  =  (m^  +  m^j .  If  this 
is  the  case,  integration  over  the  0.,  variable  transforms  the  series 
into  a  linear  ccmblnation  of  at  most  two  infinite  series  of  integrals 
of  the  type: 


00       w       /°°  /  n,+l  -k,r,,  n„-l  -k.r   /  1,  l-.+l> 


p^  (cos  e^)p.   (cose^)p^  (cos  eg3^)siii  e^de^dr^. 


-  6  - 
Passing  to  the  9^,  Sl.  variables,  we  first  expand 

^Bl%  (<^°««B1^ 

in  terms  of  r^,  Pg,  and  P^  -^  (cos  e^)  where  0  <  |iil|  <  A  <  1 

3 
by  means  of  the  expansions  given  in  Appendix  (2).  We  then  expand 

r_^      e      in  terms  of  r^,  9^,  Pg,  also  by  means  of  the 

expansions  given  in  the  same  appendix.  ISxe  result  is  a  linear 

combination  of  at  most  two  infinite  series  of  integrals  of  the 

form: 

00  % 
(3.5)     r   E   E;  M^2i_+l)pJ  (cos  a)  f   [ 

1/2  -VaI^  ,^  ,fh      J     V^\ 


X^Al 


)P^  (cos  9^)Pjr  ^   (cos  9^)P^  ^  (cos  9^) 


X  P  (cos  9^ 

^2      '^  n  J.  J 

sin  9.,  d9..  dr., 
Al   Al   Al 

where  0<  |m,|  <  /,<1,  T  =  nj^  +  lorn^  +  2  and  f^^  »  either  p^^ 

or  <i.  defined  in  Appendix  (2).  We  then  expand 
2 

(5A)     P?  (cos  a)  =  sin""  a(sin*a)  ^  (cos  a) 
^1  ^1 

M   ** 
-  Sin-  a  i:  Gi  MOi  ^  -M^Si^^^"  a)  (sin  a  ^  O). 

IsO    1       1 

This  expansion  Is  given  in  >^endix  (3)>  We  now  note  that  the 
only  possibilities  for  M  are  0,  1  or  2.  Considering  each  case  separately, 
we  can  expand 


{m.  I         l^'^l         M 
P^^  (co8  9;^)?/.   (cos  e^)P^  (cos  9^) 
131 

as  a  linear  cooibination  of  Legendre  polynomials  P.    (cos  9  .)  with 

^l""l     ^ 
the  G's  as  coefficients,  as  defined  in  Appendix  (3).  We  now  may  integrate 

over  the  9^  veiriable  using  the  orthogonality  properties  of  the  Legendre 

polynomials.  Hie  general  restilt  in  each  of  the  three  cases  M  =  0,1,  or 

2  is  that 

ft 

(3.5)     E^  „(2i2+l)p5^  (cos  a)f  P^  (cos  e^)P^  (cos  9^) 

1       Q   2  1 

kl       M 
X  py^^  (cos  e^)p^^^  (cos  e^)  Sin  e^  de^ 

=  Sin  a  2_  A    r  P.  j_  (cos  a) 
T^T   Vl  5  2  1 

vhere  t.  ,  t  take  on  only  a  finite  niimber  of  integer  values .  This  number 
is  independent  of  i^,  and  the  A.      are  rational  functions  of  i^,   vith 
the  degree  of  the  numerator  equalling  the  degree  of  the  denoninator. 
By  the  orthogonality  conditions  on  the  Legendre  polynomieJ-s,  i.,  ip 
are  now  not  independent  so  the  A.     may  be  regarded  as  a  function 
only  of  ip.  By  the  general  result  stated  above 

11m    A.  m     constant  ^  +  00,   0 

so  that  the  coefficients  A.     do  not  affect  the  convergence  of  the 

^2*^1*^5 
infinite  series  of  integrals. 

After  integration  over  9  .  the  Infinite  series  becomes  a  finite 

linear  combination  of  Infinite  series  of  Integrals  of  the  l^pe: 
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00     1 
00  (     T„+  T.   -k,r 


(3.6)    E  P,„^(=os  a)j  rl    \     1  ^n^^(.^,r^,p^)  (^pj  ^  /  „J«\  er 


Al 


vhere  the  t's  are  Rmall  fixed  integers,  and  where  the  subscript  on  i  is 

dropped.  Since  f .    »  p     or  q^    as  defined  in  Appendix  (2),  and 

-1        3      3       3 
q^    =  kp  p     plus  (higher  order  terms  in  i),(#ee  Appendix  (2)) we 

may  take  t,_^     =  P^^^  •  Writing  p.^   as  a  linear  combination  of  the 

fimctions  7j._   ^i^  "  f-  +  1,  we  can  expand (5. 6)  as  a  linear  combination 

of  at  most  two  series  of  the  following  type: 

(3.7)     E  Pi+^^(cos  <=^)(2lT2^  S^ 

OD      3    »- 
®  1      r  '^2'*"  2  "^1* 

-  \-   ^i+T^(^°«  ^^  (2i-H2T+l)  J  ^       ^     ''i+Ti^^^2'^'^B^ 


'O 


X  {  Pp^  /  Pc  >  )  *^     where  x  =  r 


Al' 

We  win  now  proceed  to  obtain  the  rate  of  convergence  of  the  8erie8(5-7) 
and  the  asymptotic  behavior  of  the  terms  of  this  series  for  large  i. 

The  general  procedure  is  to  split  the  region  of  integration  into 
subregions  where  the  integrand  is  anedytic  and  to  consider  each  term 
of  the  series  as  a  sum  of  integrals  over  each  subregion.  Then  each  of 
the  summands  is  cmalyzed  sepeirately  and  in  this  way  the  overall  behavior 
of  the  series  is  obtained.  Each  suimwnd  is  analyzed  by  first  noting 
that  the  integrands  are  non-negative.  Ihen  if  the  integrands  are  re- 
placed by  non-negative  functions  which  are  never  less  than  the  integrand 
and  the  integration  carried  out,  a  majorant  series  for  the  original 
series  will  be  obtained.  Similarly,  by  using  functions  which  are  never 
greater  than  the  integrand,  a  minorant  series  for  the  original  series 
will  be  obtained.   If  the  majorant  and  minorant  series  have  the  same 
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asymptotic  behavior  and  the  same  rate  of  convergence,  we  may  then  con- 
clude that  these  are  the  asymptotic  behavior  and  convergence  rate  of 
the  original  series.  The  functions  used  to  replace  the  integrand,  to- 
gether with  other  inequalities  used  are  derived  in  Appendix  (k) .     To 
carry  out  this  procedure  we  must  consider  three  cases  Pp  <  p_,  p_  >  Pg, 
and  pp  =  Pg. 

Case  1.  p^  <  pgt 

If  we  ignore  the  coefficient  to  the  left  of  the  integral  sign 

for  the  moment,  S.  can  be  written  as  a  sum  of  three  integrals: 

12    5 
^i  =  ^i-^^i^^i 


where 


a    -i-i  r ""  ^2  *  2  -^  ^  -V 


(5.8)     ^1  =  ^c    J  ^         ®     ''i+T  (^2'^''^B^  ^ 


Pb  1 

2    i  /    ^2  *  2  "  ^  " V 

(5.9)     S^  -  Pc  J   X  ^  ^     e  -^  ^+Tj^(V^>Pb^  ^ 

Pn 


5    w'^  T  +  ^  -  i  -k  X 
(5.10)    S^  -  pj  J  X  ^   ^    ^  ^i+r    (V^'Pfi^  ^• 

^B 


We  may  now  use  the  inequalities  of  Appendix  (k)   to  deduce  the 
following  inequsilities : 


"Vb  -(V^2^Pc  -\-   I  ^2^\ 
1    ®      ®          ^B      ^C        "^C  i 
^^•^^    ^i^     (2i+T2+Ti^+5)(2i+2Tj^+l) (^) 
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1   .^  ^B  PQ  Ai 


^^•^^^  ^i  ^       (21+T,+T,+3)(2i+2T^+l)  (^) 


-k«P^  -T.  -  -^ 


3  /^Cxi    f         ^a'^i**^     -(lE3^+k2)x 


'^B  "'k'  2  /b 

^^•^5^  ^1^     (2i+2Ti^+l)  ^^M        ""  *  ^ 

_      -  1 


^  V^l^^B  ^2*  2 


5  ^  ^B  /C  i 


^^•^^^  ^i  ^  (2i+Ti^-T2-l)(2i+2Ti^+l)    ^^^ 


Case  2,  Pg  <  p^: 


If  we  agala  ignore  the  coefficient  to  the  left  of  the  integral  S. 
can  again  be  written  as  a  sum  of  three  integrals: 


12         5 


where  nofw 


p. 


B  3 


(3.16)         s\    -     p'^"^     {       x^         ^  ^     ^     ^1+T  (V^^'^B^  ^ 

0  ^ 


p  -1.1        A      T   +i+  I      -k^X 

(3.17)        sj    =    p^^  ^   J      X  2        2  ^  ^    ^+t^(V^>Pb^  ^ 

^B 


T-,+  ;;  -1     -k,x 


(3.18)        s5    =    pj  X  2    2        e    ^    r^^    (^x^pg)  dx 


^C 


The  following  inequeilities  may  readily  he  deduced  from  Appendix  (k) , 


11 


.(k^+2k2)P3  T2+  I  .1 


Pn.i 


^^2-^1 1*^8  ^2-^  2  -1 
1     ®        Pg    Pg     Pg  i 

^^•^°^    ^1  -  {2i+2T^+l)(2i+T^+T2+5)  ^^^ 

-k  p  T.  +  ^  pp 

(5.21)    S?  >  ,^,^^,1,,       ^   (-2)     X  2  ^   e   ^  2   ^ 


'i  -  (2i+2Tj^+l)       ^p 


C 


•^B 


X       e    dx 


^2^B  "^U*  2  -1         ^C 
_     e  ^  ^^^  2p/   p_  i  ("^     T.-T,.+l  -k,x 

(3.^)   sf  <  (2i4.i)    (;^)   ^ 

Pb 

5     e     e     Pg    p^  Pg  1 


i  - 


(2i+2Tj^+l)(2i+Tj^-T2-l)      ^p^ 


Figure  2  summarizes  the  asymptotic  behavior  of  the  S^,  J  =  1»2,5 
in  both  the  preceding  cases  deduclble  from  the  above  Inequalities. 
Returning  to  the  original  series  (5.7),  and  now  taking  into  account  the 

coefficients  heretofore  neglected,  we  see  that  the  series  converge 

C   ^<  i 
absolutely  as   )     ( — )  for  large  1.  "Hiough  we  only  obtained  a 

P'       P> 
majorant  series  for  S^,  this  is  sxifficient  for  us  to  determine  the 

asymptotic  behavior  of  the  series(3.7). 

Case  3»  Pg  »  P(,: 

Here  S.  =  0  and  S  <=  S  +  S^,  where  S  ,  S^  can  easily  be  seen 
to  behave  asymptotically  as  (T/i  ,   C-^/i  respectively.  Hie   absolute 
asymptotic  behavior  of  the  series  (5'7)  in  this  case  is  „  i^   for 
large  1. 


.1  i:, "-<)'■ 
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1    i  ^pJ 


s5  <  ^(^)' 


p^  =  min(pg,p^) 


p^  =  max(pg,p^) 


12  3 
C  ,C  ,C     are  fixed  constants. 

FIGURE  2 


1|..       Two  center  exchange  Integrals 
Biese  are  of  the  form: 


(4.1)  j    (  ^[^\lhi^\l)r-^^[^h2H^^^\2)  dv^dv 


2  1 


00    00  n     n     2k     2it 


^'-^1    -^3''a2 
e 


X 


^C^C^CC     r     ^     D        (     j     j     j     J        J 

V^     -VaI  V^     -^2^B1,    ^    /  /+1n3' 
X  r^^;'  e''^'«2  p]^'(cos  e^)p|^'(cos  e3,)^(coB  9^) 

X  ^(C08  e^)P^"^  (cos  Q^)^j^  (cos  ©02)^^"  ®A1^^  ®A2 
\^^ja^  \%l^'''>^^^Al'^P2^  ^^^A2>  ^(i^B2^'^V<^A2'^®Al^®A2 


dr  .dr 

Al      A2 


15 


^  \ 


Beginning  with  the  0  ,  0^  variables  we  first  expand  ^  ^  and 
each  as  a  linear  ccanbination  of  sines  and  cosines.  The   integral 


vill  vanish  iinless  at  least  one  sine  or  cosine  fiinction  in  the  expansion 
of  3   ^  is  the  same  as  a  sine  or  cosine  function  in  the  expansion  of 
i^  C.  .  In  this  latter  case,  integration  over  0^,  0^  transforms (4 . 1 ) 
into  a  linear  combination  of  at  most  two  series  of  the  type: 


00  00  rt   : 


-  -  «  -«  n,.l  -k,r^  n,.l  -^31 
^Al   ^      ^Bl   ^ 


0  '0  0 
X  (^A<  /  ^A>  ^  ^A2   ^      ^B2   ^ 

X  pP-  (cos  e^)P^  (cos  eg^)p"(co8  ©^)p^(co8  e^) 

|m_  I  [ffli  I 

X    Tp'  (cos  e^)PX       (cos  eg2)sin  9^  sin  9^ 

X     de^  de^  dr^  dr^. 

Passing  to  the  variables  9^,  9^  we  expand  r^^  P^   (cos  9g^)  and 

TgJ  P£   (cos  9^2)  in  terms  of  r^,  p^,   P^  ^  (cos  9^)  0  <  |m^|  <  l^  <  1 

and  T^,   pg,  p|°^'(cos  9^)  0  <  |m^|  <  /^<1,   respectively,  by  the 

formulas  of  Appendix  (2).  Next  we  expand 

»2-  ^2-^   -^2^B1 
^Bl 

eind 

n^-  \-l       -kj^rgg 
^B2 

in  terms  of  r^,  P^(cos  9^),  pg,  and  r^,  P^  (cos  9^),  pg,  respectively. 

(See  same  appendix)  Ibe  result  of  these  expansions  tran3forms(J<-.2)  into 

a  linear  cooibination  of  at  most  four  series  of  the  form: 


Ik  - 


00   00   A   It 


OD   00    00  f         (  I       ( 

(^•3)     E  E   n  D^j/2i  +l)(2i  +1) 

X^A2    *      %(V^A2'^b\(<^°«  V^/^(=°^  V 

X  p^^  (cos  e^)p^(cos  e^)p^(co8  ©^)p^^  (cos  e^) 

X  P^^  (cos  9^)  P^^(cos  e^)8in  e^sl»  e^de^dG^dr^dr^ 

where  t,,  t^  s  1,2,  or  5,  and  f .  =  either  p.  or  q.   defined  in  Appendix 
(2). 

Now  M  may  take  on  only  one  of  three  possihle  values:  0,  1,  or  2. 

M  KI 

Considering  each  case  separately  we  may  expand  R  (cos  0„,)P/   (cos  6.,) 

KI  M  KI  l"»6l         ^ 

X  P^  -^  (cos  e^)  and  P^(co8  9^)?^^   (cos  ©^)P£   (cos  6^)  each  as 
a  linear  conibination  of  Legendre  polynomials  with  the  O's  as  coefficients 
as  defined  in  Appendix  (?)•  Ihen  the  integration  may  be  carried  out  over 
®aT  ®a?  ^y  *^®  orthogonality  relations  between  the  Legendre  polynomials, 

Bie  general  result  in  all  three  cases  Ma  0,1,  or  2  is  that 

n  *  11 

{k.k)  D^j/2i^+l)(2i2+l)  I  J  P^  (cos  e^)P£^  (cos  e^) 

KI  _M  KI 

X  P£  ^  (cos  e^)p"(co8  e^)p^  (cos  e^)p^  ^  (cos  e^) 

5  2         3 

KI  M 

X  p^^  (cos  e^)p"(co8  e^)sin  e^sin  e^de^de^ 

6 
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where  A /7     is  a  rational  fiinction  of  /..  Hie  degree  of  the  numerator 

is  the  same  as  the  degjree  of  the  dencmiinator,  and  the  upper  limits  on 

the  svmimations  over  t,,t.  are  Independent  ot  iL  .     By  the  orthogonality 

properties  of  the  Legendre  polynomials,  i^,  i^  are  no  longer  Independent 

but  become  functions  of  /C  •  Since 

lim   A  /;     ■  constant  /  -f  oo ,  0 
/->oo  ^V»^ 

by  the  remarks  above,  the  coefficients  A /?     cannot  affect  the  basic 

convergence  of  the  series.  So  integration  over  0^,  Q^   reduces  the 

series  to  a  linear  combination  of  series  of  the  following  type: 
OD  00      1   , 


0  0 


X  (^A<  /  ^>  )  ^A2     ^      ^/+  T  (H'^A2'Pb^ 


X  dr  ,  dr  „ 
Al   A2 


As  in  Part  3,  we  may  take  t ^    :  p.  (Appendix  (2))  and  write  v />    as  a 
linear  combination  of  7/f.       (Appendix  (2)).  lhen(4.5)reduces  to  a 

linear  ccmibination  of  at  most  four  series  of  the  type: 

•  QD  -  _.^  1  ^. 
00         -  f     f       T^+  p  -k^r^ 

^""'^^  ^  (2/+2T,+l)(2^+2T.+l)  J   J   ^Al    ^ 

^^+T5<V^A1'Pb^^^A<  /  ^A>  ^^A2    ^ 

74Tg(H'^A2'PB^  ^Al'^A2 
where  t^  r  t,  +  1,  t^  ■  Tk  +  1. 
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noting  the  definition  of  (r^  /  r^>)  we  write  (4.6)  as  a  sum  of 
two  series,  (h.l)   and  (k.Q)   where  (4.7),  (4.8)  are  given  by: 

•Al 


00  ,         r^  Fp2    {'+t-,^+  I  -kj^r^ 


^^''^^  ^(2i+2r5+l)(2/+2T^+l)J  J    '"Al 

X  7^+-  (V^M'^B^  ^A2 


X  ^/.Hg^H'^AS'^B^  ^Pl^t^ 


e 


(1* 


a>  «>   i+T,+  i  .k,r^ 


■®^     ^  (2A+2t,+1)(2/+2t,+1)J   J   ""aI 
^  0   r^ 

>^+T2+  2  ■^3^A2 
X  '£+T^(*^2'^Al'^B^  ^A2 


X  7/4^^(^>^A2''^b)  ^A1*^A2 


If  in  (4.8)  we  Interchange  the  order  of  integration, 
00  00     00  r., 

0   r^     0   0 

we  transform  (4.8)  into  a  series  of  the  type  (4.7)  so  we  need  only 
consider  series  of  type 


00   X 


00  -^  00  ^  f  f 

^''•^^     ^  (2i+2T^+l)(2/+2T^+l)  \  '  f   (2/+2T^+l)(2/+2T^+l)  j   j 

/+T  +  I  -k^y  -  ji+T  +  I  -k  X 
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vhere  y  =  r^,  x  -  r^. 

Ignoring  the  coefficients  to  the  left  of  the  integral  signs 

12    5 
momentarily,  ve  vfrite  S />  as  a  sum  of  three  integrals  Sy  «  S2  +  S^  +  3/ 

where  the  integrand  of  each  Si ,   J  =  1,2,5  is  analytic  over  the  region 

of  integration. 

P 


(4.10)  4     "    j       J     y         ^         e     ^ 


0      "0  ^ 

-  I+T2+  I     -k,x 

X    X  e     -^     7^^^  (kj^,x,Pg)dydx 

GO      "^B      ^  5       , 

(4.11)         S|     -    j      j       y         ^     2  e     ^     WV^'^b) 
Pg     0 

-  Z+T  +  I     -k  X 

XX  '^     ^  e     ^     7^_j^  (ki^,x,pg)dydx 

°°    -/4a_+  I     -k  X 

Pb 
P 


0 


x]   y    "  ^e  -  ^^^^cv^^pb)^- 


/°°  /    /.-«,+  I     -k  y 
(lf.l2)         SJ     -     j     j     y  '^  e     ^     7^^    (k2,y,PB) 


Pb^b 


X  X  e     ^     7.^   (k|^,x,PB)dyiix 


°"   °°      >^-^2*l     -V 


J    /    X  -     -  e     -^     7^^^(k^,x,P3)dy 


^B^ 

-k^y 
X     «  ''/»j,   (ko>y>pTi)clxdy. 


'£+r^''2''''^B' 
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Proceeding  as  before,  by  the  Inequalities  of  Appendix  (k) ,   we 
can  readily  deduce  the  fo]J.owing: 

1        Pfi      ® 
^^'^^^         ^i    -     (2/+  ^5+l)(2Z+&r^l)(2X+T3^+T5+3)(2£+ri-Hr5+T^5) 


Pt,      e 


1  B 


'^l+'fg"^!  (kj^-k^+(k2-k^)+)pg 


2         *^B      " 
^^•^^^    ^l     -     (2X.+2T^+l)(2/+T^+T5+5)(2/+2T^l)(2X-T2-W6-l) 


,.  .       B    ^ 

'/-  -  (2/,+2t  +l)(2i+2T^l)(2/-T2+Tg-l)(2/-T2-T^+T  +T^1) 


(i*.l6)    S^  -  ^ 


Hence  S^  can  be  both  majorized  and  mlnorized  by  terms  of  the  form  Cr/^ 
(:  o  constant  for  large  £, .  Ihus  S^  behaves  asymptotically  as  Cr/^ 
for  large  /  .  Since  both  Sf,  Sj  are  majorized  by  expressions  of  the 
form  Cr//.  ,  C^/ a  ,  where  C^,  C^  are  constants  it  follows  that 
Si  »  S^  +  S^  +  S^  must  behave  asymptotically  as  C//  for  large  ji.  , 
C  »  constant.  Therefore  it  follows  that  the  series  (^.9)  must  behave 
as  c/£.  for  large  jt  if  we  take  into  account  the  coefficients  to  the 
left  of  the  integral  signs. 


-  19  - 

5.   Three  center  hybrid-exchange  Integrals 
Ihese  are  of  the  form: 

(5.1)  I    j     t^^^(l)Mr^^^(l)r-^tp)(2)4''^(2)dv^dv2 

2     1 

C»    QD     )l      Jt      2«      2jt 

=  c^c  c  c^  f;   ;^     °/m  f    f    [    (    f      f 

^2  5  4^^Q^0  i     i     I    i    i       i 

„,.!  -V       n,.l  -k^r^  ,V\-^3^A2 

Al  Bl     ^  ^^A<  '      A>    '      A2     ® 

n.-l     -k.r         luL  I  |nu| 

(dl  I  |ini  I 

X  ^(cos  ©^)PjP  (cos  e^)Pt     (cos  ep2)sin  e^sln  9^ 


i,(^Al)  i^^i^Bl)-^  "^(^A1-^A2)  i^(^A2)  i,(S^C2) 


X  d0^  d0^  de^  de^  dr^  dr^ 


We  first  expand 

|n»6l 
in  tems  of  r^,   p  ,   and  P£      (cos  9^)  ^  (^'a?^  ^^  *^®  expansions  given 

6  6 

in  Appendix  (2).  nien 

L       |mJ 
^Bl  V^^°'Sl^ 

|^'  ,   .    . 

is  expanded  in  terms  of  r^,  p_,  and  P^  -^  (cos  9^).  Here  0  <  |ai^|  <  £^  <  1  , 

.    .    ^   «  .        ^  VV^  H^Bl   VV^  '^02^  ,^ 
i  =  5  or  6.  Next  we  expand  r_^       e     ,  r-^      e      by  the 

expansions  given  in  Appendix  (2).  The  series  (5-1)  then  becomes  a 

linear  combination  of  at  most  six  series  of  the  form: 
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(5.2) 


00        X       00  CD  2  J_ 

r    r    E      L      E       I^£i21  +!)£  P^  (cos  a) 

jT^O  n^O  i^=0  1^=0  Jg«0    '^  2"^  2     ^2 


oooort«2jt2K,.l       ,^ 

0     0     0     0  0       0 

^  f    i     /     t+1^     '^2*  2     "S^A2   .     /.  . 

kl  M 

X  T^     (cos  e^)P/     (<=os  e^)p^(cos  e^)p^  (cos  9^) 

X  Pjj^  ^  (cos  e^)p^  "^  (cos  e^)iJ'(cos  e^)p^  (cos  e^) 
36  2 

X  sin  e^sin  e^  ^(0Ai)  \^^i^^  <=°«  °(<^A1-^A2^ 
X  ^(i^A2)  ^(^A2^   *^°«  J2'^A2  ^Vl^A2  '^V^®A2 
X'^A1^A2 


where  t  ,  t  =  1,2  or  3,  f^  =  p^  or  q^^,  as  before. 

We  expand  <B   ^  and  ^^  ^  cos  J-  0   each  as  a  linear 
conibination  of  sines  and  cosines,  liie   Integreil  vlll  vanish  unless 
a  sine  or  cosine  fimctlon  in  the  expemslon  of  \^  \^    equals  a  sine 
or  cosine  function  in  the  expansion  of  5^  ^  cos  J^  0^2'  ^  *^® 
latter  case,  integration  over  0^,  <f)^   transforms  (5«2)  into  a  linear 
conibination  of  at  most  two  series  of  the  form: 
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(5.5)     n  r  zi  i^^M^^v^^^i  i^i  ^^°^^^ 

£=M  i^r=0  ±2=  J   ^M   -L    ■i-2'J  -"2 


0  0  0  0 


^B'^"A<  '  ^A> 


X^A2  ^  ^   V^'^'"a2'^C^V^'°'V%^^°^V 

2  15 

X  if(co8  e^)p  (cos  e^)p^"5  (^^3  e^)p  "6  ^^^^  ^^j 

15  o 

X  p^(co8  e^)pJ^(co8  e^)8in  G^sm  e^  de^  de^  ^ai"^a2- 

As  before,  M  =  0,1,  or  2  only.  Then  J  «  0,1,2,5,  or  k   only.  Similarly, 

we  expand  P'J  (cos  a)  =  8in~''a( sin*'  a  P^  (cos  a))  as  a  linear  comibination 

ig  ^2 

of   Legendre  polynomials  vith  the  G's  as  coefficients  (see  Appendix  (3)) 


By  the  expansions  in  the  same  appendix,  P£   (cos  0^)  T£   -^  (cos  e^)ly(cos  0^) 


may  be  written  as  a  linear  combination  of  Legendre  polyncmilals  with  the 

G's  as  coefficients.  Next,  we  expand  P^   (cos  ©^)P/   (cos  G^) 

X  pJ^(co8  G  ,JP^  (cos  G„„)  as  a  linear  combination  of  products  of  two 
Jt^     ii2  ip     A2 

Legendre  functions  with  G's  as  coefficients.  The  analysis  of  all 
possible  cases  arising  from  M  =  0,1,  or  2,  and  J  =  0,1,2,3  or  h   show  that 

i^.k)  8in"'^a( sin-'a  pj  (cos  a))D£j^2i^+l)E^  . 

xjj    P^"^  (cos  G^)p/5  (eo8  G^)pf(cos  G^) 

0  0  ^ 

|m,  I         |m^| 
X  P^  (cos  ©^)P£   (cos  G^)P£   (cos  G^) 

X  P^(co8  0^)P^  (cos  G^)sln  G^sin  G^  ^G^^G^ 
-  sln'^a  E       A  T..T  T  T  ^/+T  ^°°^  *^ 
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where  the  A  /;  are  rational  functions  of  i.   the  numerator  having 

the  same  degree  as  the  dencanlnator .  The  upper  limits  on  the  summations 

over  the  t's  are  independent  of  /-   .     By  the  orthogonality  properties 

of  the  Legendre  functions  the  i^ .  i-  are  no  longer  independent  but  are 

functions  of  L  .     Since  by  the  above  remarks  lim    A.         =  constant  /  0 

^  ->  00  ^'^8'^9  lo"^? 
i  +  CD ,   the  coefficients  A/,        do  not  affect  the  convergence  rate 

•^'^8^9  lo'^T 
of  the  basic  series  of  integrals. 

Integration  over  the  0^.,  0.-  variables  reduces  (5«5)  to  a  linear 

combination  of  series  of  the  following  type: 

CO  00     1   _. 
(5.5)     r  P^   (COS  a)|  f  r^  "  e"  ^  "^^/.^  ( V^Al'^B^ 

X  (-f<  /-A'>')-Ii*  '  «'''^''^^^^JV-A2^Pc)  '^A1^A2- 

As  with  two  center  exchange  integrals  we  note  the  only  possibilities  for 

the  f  functions  and  write  (5-5)  as  a  linear  combination  of  four  series 

of  the  type: 

°°         1  [     (        1  2  ~  I'^Al 

^5-6)     ^  (2i-^2T,fl)(2/^2T.4l)  ^/^T„^"'^"  ^M  j   ^Al    ^ 
-^^         ^  '       0  0 

><^£+T5^V^Al'^B^(^f</^''>)^A2  ^^  ^  ^^i+T^(H'^A2'Pc^ 

X  dr^  dr^   where   T^  =  t^  +  1,   T^  =  t,^  +  1. 

Noting  the  definition  of  (rZ;  /  ^a-s.)  "^  write  (5«6)  as  a  sum  of 
two  series,  (5.7)>  (5.8)  where  (5-7),  (5-8)  are  given  by 


(5.7) 


(5.8) 
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00  ^A2 


Pj?4-.    (cos  a) 

^^1    

(2X+2T  +1)(2£+2T.+1)        11  ""aI 

0     0 


C  +T  +  ^     -k  r 

^      1     2       ^1  Al 


><  ^/+T  (V^A1'Pb^^A2       ^     ^  e     5  ^7^^     (ki^,r^,Pc)<ir^<ir^ 
5  o 


P^^^^(co8  a) 


(2/+2t^+1)(2X+2t^+1) 


00    00  «  1         , 

.         -^+V2     -VaI 


0     r 


A2 


jI+T  +  ^     -k  r 
^  ^/+T  (V^A1'^B^^A2  ^     ^  e     5  A2  7^^^^(k^,r^,p^)dr^dr^. 


Interchanging  the  order  of  integration  in  (5-8), 

oo  .00      00  r. 


0   r 


A 


Al 


A2 


0   0 


transforms  the  series  (5.8)  into  a  series  of  type  (5.7),  so  we  need 
consider  only  series  of  type  (5.7)*  namely. 


(5.9)      n 


OD       i+T 


T.        (cos  a) 


j-     (2X+2t^+1)(2^+2t,^+1)  ^X 


^  P.^     (cos  a)  00  X       ^^     ^  3       , 

00  /+Tt  f     (        £+T3^+  I     -k^y 


y-      -  ^1   ^  f  ( 

"  ^  (2,£+2tj+1)(2/+2ti^+1)  I    )     ^ 


><^^+T5^Vy''^B^^ 


0  0 


X  dy  dx    where  x  =  r^  ,     y  =  r 


AL  • 


We  now  must  consider  three  cases  p_  <  p_,  p^  >  p„,  and  p_  t=  p 

D  \j         D  \j  a  \, 


Case  I,  pg  <  0(, 


Ignoring  the  coefficient  to  the  left  of  the  integral  signs,  we 
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write  S£  =  S^  +  S^  +  s|  +  S^  +  S2.  +  S^  where  in  each  s4  the  integrand 
is  analytic  over  the  region  of  integration. 


(5.10)    5=1   I  y      e   7^+T  (k2,y,Pg)x 


'0   '0  ^ 


X  7^^  (kj^^x,Pp)dydx 


^C  -Z+r^+  I     -Kx 


(5.11)    S^  =   I   X        e  ^  7^^  (kj^,x,p^)dx 


^B 


X  1   y      e    7^^^  (k2,y,P3)dy 


0  5 


,^C  /  i+T  +  I  -k  y 
(5.12)    SJ  -  J   J  y       e  -^  7^^^  (k2,y,PB) 


f^B  ^B  5 

-I+T2+  I  -k,x 
X  X        e    7£^  (k|^,x,pp)(iydx 

6 


^2,     r°°  "/+T2+  I  -kjX 


(5.15)    ^  =  I  X    ^  ^  e  ^  7^^  (kj^,x,p^)dx 


^C 


6 


,^B/+T+|  -ky 


•0  5 


"^   -i+T.+  i   -k. 


(5.1^)    ^5  =  J   X    ^  ^  e  5  ^_^^  (ki^,x,P(,)dx 

/^C  j^+T  +  I  -k  y 
X  j   y       e  -^  7^+T  ^^2'^'^^^^ 
^B  ^ 


-  25 


a>  00./+T  +i     -V  X 


(5.15)  ^      "      1     J   "^  ^         '>+T.^H''^'Pc^ 


P^y 


X  y  e         7£^T  (k2.y.pB)<ix^- 


From  the  Inequalities  of  Appendix  (4)  we  can  readily  deduce  the 
following: 

-(k    +2k    +k    +kj^)p3    -k^^P^   I   +T^+T^   +T       -T^-    i 

1  ®  e  Pg  p(,  Pg^ 

(5.16)  S^>      (2^+2T^+l)(2^+2Tg+l)(2/+T^+T^f5)(2^+T^+T2+T^+T5+5)    ^^^ 

(k^-k  +(k2-k^)^)^P3  I  -Ki+T^-Hr     -  Tg-  I 

.  .  1  ^  Pb ^C^ ,%^ 

(5-17)  ^  <  (2£+2T  +l)(2X+2Tg+l)(2X+T^+T^+3)(2;e+T^+T2+T^+Tg+5)^Pc^ 

-(k    +2k    )p  -kj^p^    T   +  I   -Tg-   I  ^      .    ,PC    ^   ^^   ^-        ,,     ^,      ^ 

2  e  e  Pg^       P(^  "^B  W       ^2^6       ■(^5'''H^ 
(5-18)          S^>     (2/+2t6+i)(2Z+2t  +1)(2/+T^+T  +3)   ^T^  )     ""  ^  '^ 

Pb 
(k2-k^)PB  Ti+  1-^6-1  „  ;e  ,^C      '^  ^,   ^-     .    ^ 


^5-19)  S^  <  (2^+2Tg+l)(2/+2T^+l)(2X+T3^+T^+5)   ^Pc^  j     "" 


(_£)    1     X  -^     ^     e     -     -^     dx 

Pt 


^B 


Vb    "^^Pc  '^5"'  2  "  "^6-  2  ^ 

e  e  p„         p„  Pt^-*^ 


^  B  C  B 

^5-20)  S^  >     (2^+2t5+1)(2£+2t^1)  ^^^ 


/"C    f       T,-T.+l   -(k  +k    )y  T  +Tg+1  -(k  +k.  )x 
xj       Jy-^^e      ■^'^x'^^e      ^    ^      dydx 


Pb   Pb 


T^+  o^  T^-  t: 


^5"  2  ^6    2  . 

3  ^B         ^C  ^B-^ 

(5-21)  S^<     (2£+2T  +1)(2/+2V1)   ^r^ 


Pc   ^     *'i-''=+l     -(k,+k  )y  Tp+^6+1  -(k.+ki^)x 
X  y^       e       ■^'^x'^^e       ^     ^     dydx 


Pb    Pb 
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(5.22)  ^/  ^  (2/+2T^+l)(2/.+2Tg+l)(2/+T^-T2-l)(2/-K;^+T^+3)    ^0^,^ 


Vb  ^H-^3^^C  ^'f  2  ^2-'  2 

e         e  ^       P_         p„  p 


B 


-     5 
y  e         dy 


/ 


5  B  C  B 

(5.25)  S^  <   (2/+2T  +l)(2X+2Tg+l)(2X+Tg-T2+l)      ^H^^ 

'Pt 

Vb  (H-^1-^2)Pc  •^5*  i  ^l-^2-''5"'  2  / 

gee  Pg         p^  ^  Pg  ^ 

(5.24)  S^  <  (2^+2T^+i)(2/+2Tg+l)(2;(.-T2+Vl)(2^+Tg-T2+T^-T^-3)    ^^^ 

We  sxamnarize  the  asymptotic  behavior  of  the  S^  for  large    X    as  a  resvilt 
of  these  inequsilities  in  Figure  5.     Fran  this  we  deduce  that  the  asymptotic 

c     Pb  -^ 

behavior  for  large  JL     of  S/»  is  as  —  ( — )  C  =  constant.  Hence  the 

X.  £^      Pp 

absolute  asymptotic  behavior  for  large  X     of  the  terms  of  series  (5.9)> 
taking  into  account  the  coefficients  to  the  left  of  the  integral  signs 

c   Pb/ 

is  as  -T-pr   (— ) 


"972 


P 


C 


Case  2;  Pg  >  p^. 

We  again  vrrite  S^    =  S£  +  s|  +  S£  +  S^  +  S^  +  S^  where  now 

(5.25)  S^     =     j       Jy  e-"  7^_^^  (k2,y,PB) 


0      '0  ^ 

/+T2+  I     -kjX 


X  X  '^     ^  e     ^  7.^.^  (k^^,x,pp)dydx 


Pt 


(5.26) 


S|    =    J     X  ^     ^e     ^  7^       (k^,x,p^,)dx 


Pc 


/^J^+T+I     -k^y 
X    J     y       "•         e         7^^     (k2,y,PB)dy 

0  -^ 
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y  — > 


The  C*^  are  constants. 

The  double  lines  bound  the  region  of  integration; 

0  <  X  <  +  00, 
0  <  y  <  X 


FIGURE  3. 


(5.27)        s5    =    j      j 


27  - 


^c    ^C 


XX  e     ^  7£^  (kj^,x,Pj,)dydx 


<»  £  +T^+  i     -k.x 


(5.28)  Sj     =     J        X  e     ^  7^^  (kj^,x,p^)dx 


/  "  £+T  +  §     -k  y 


^  f^    -  £-H-  +  I     -k  X 

(5.29)         s|    =     1     X  ^  e     -"  7^^  (ki^,x,p^,)dx 


Pb 


.^B 


(5.50)  S^    =     I       J     X  '^     2  e     ^  ^;+Tg(^'^'Pc) 


^B  y 


;^+T3^+  I     -kj^y 


X     y  e  7.^  (k2,y,PB)dxdy. 

A     5 

From  the  Inequalities  of  Appendix  {k)  we  can  readily  deduce  the 

following : 

-^2Pb  -(ki+2k2+k  +2k^)pp  -T  -  I  T^+T2-»-Tg+  I 

1  ®  ®  Pg  P(,  PQ  ^ 

(5.51)     ^1  ^  (2X+2T^+l)(2/+2Tg+l)(2/+T^+T^+5)(2/+T^+T2+T^+T^5)  ^P^^ 


(5.52) 


si< 
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^c/ 


(r^) 


k   -   (2^+2T    +l)(2/+2Tg+l)(2£+T^+T^+3)(2X+T^+T2+T^+T^+5)       Pg 


•(  V^5+H)p'B  -( VVH^^c  -'^5'  2 1  +V^2-*-^^ 


2  '^B  "C  /[X 

(5.55)  ^X-  (2X+2T  +1)(2/+T^+T  +5)(2X+  2T^l)(2iL+Tg-T2-l)      ^Pg 


^c/ 


X 


Pp    2£+T,-T    -1 


(5.5M 


((k2-k^)^-k    )p^-T    -||+T^+VT. 


s!"  < 


^B 


^c/ 


(7^) 


'X  -  (2£+2T^+l)(2X+T^+T^+3)(2^+2T^+l)(2£.+Tg-T2-l)    ^Pg 


X 


(5.55) 


(5.56) 


"Vb  -(ki+k2+2k5+ki^)p^  Tg+  I  -T^-  I 


sl> 


^B 


(-) 


(2X+2T^+l)(2£+2Tg+l)(2X+T^+T^+5)  ^B 


^B  "C 


T^+T^+T^-Tg+i. 


(2£.+Ty-Tg-l) 


p       2L-T    +Tg-1 


S-;  < 


(k   -k    )^Pg  (ki^-k2+(k2-k^)^)p^  T^+  I  -T   -  I 

e  e  p^     Pg   ^  ^p^^^ 


'/-         (2X+2T  +l)(2/+2Tg+l)(2£.+T^+T^+5) 


■'B 


X 


^B  -^C 


T^-K^-^^-Tg+l. 


T2Z^VV^ 


Pc    2£ -Tg+T^-l    ^ 

^B 
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(5.57)    Sj  <  (2X+2T^+l)(2£+T^-T2-l)(2/+2T^+l)(2/+T^+T^+3)   ^^^ 

(kj^-(k2-k^).)p^  (V^l^+V^3^  ^B  V  2  V'^2-'^5"V^ 
c   e  e  P^,    Pg  p^  X 

(5-.58)    S^  <  (2^+2T^+l)(2^+T^+3)(2/.+2Tg+l)(2/+Tg-T2-l)         ^~^ 

(k2-k^-k^)p3  kj^p^  ^i+Tg-V  I  ^6-^  I 
g       e  e    pg         p^  ^Pj^  ^ 

(5.59)    Sj^  <   (2X+2T^+l)(2/.+2T^+l)(2/-T2+Tg-l)(2/+T^+T^-T^-T2-5)  ^P^^   • 


From  these  Inequalities  we  can  conclude  the  following  asymptotic 
for  large  X-  summarized  in  Figure  h . 

From  this  we  conclude  that  the  asymptotic  behavior  for  large  ^^ 

C   ^C 
of  S,  is  as  —  ( — )  C  =  constant.  Hence  the  absolute  asymptotic 

A     X^  Pg 

behavior  for  large  ^  of  the  terms  of  the  series  (5-9)  is  as 
P,^ 


— ■■■■  /  ( — )   taking  into  accotmt  the  coefficients  to  the  left  of 
^11/2  Pg 

the  ingegral  signs. 


Case  3:     Pg  =  P^, 

2    5    5 
Here  S.  =  S^  ■  S^  *  0;  and  hence  the  absolute  asymptotic  behavior 


of  the  terms  of  the  series  (5.9)  can  be  shown  to  be  as 


29a 


^B 

^C 

//                 i         B 

/ 

c      p5    Pp  / 

A 
X 

'^     /        "^B 

^B 

f 

y     > 

The  C*^  are  constants. 

The  double  lines  bound  the  region  of  integration; 

0  <  X  <  +  oo; 
0  <  y  <  X, 


FIGURE  k. 
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6.   Three  center  coulomb- exchange  Integrals 
Tliese  are  of  the  form: 

(6.1)     j  j       4^)(lK^2)(l)  r"^  *p^(2)4''^2)  dv^dv^. 
2  1 

The  analysis  of  these  integrals  is  identical  to  that  of  the  hybrid  ex- 
change integrals.  Ihe  results  of  the  analysis  are  the  same  as  those 
for  hybrid- exchange  integrals  except  in  the  case  in  which  \ifg  (l)  =  \|fg  (l) 
-  '!/^  r.  r.  -uCl)  (see  Appendix  (l)).  Here  the  reduction  of  the  integrals  to 
the  series  (5-2)  is  identical  to  that  of  the  hybrid-exchange  integrals, 
but  in  series  (5.5)  f/?    no  longer  eqimls  p.^  or  q.  ^  ,   but  equals 
l;      .         (see  Appendix  (2)).  The  expansion  of  ^  .^^     in  terms  of 
y.         gives  rise  to  linear  combinations  of  7's  of  the  form: 

^^•^^  T2Z+2t5+1)(2£+2t  +1)  ^£+T5  • 

So  an  extra  factor  of  ( ^  ,  ,^ — ttt  ioto  the  series  of  (5-9)  is  introduced. 

A  summary  of  the  asymptotic  behavior  of  the  series  (5-9)  in  this  case 
is  given  here: 
1.   Pb  <  Pc 

,11/2 


The  terms  of  series  (5-9)  behave  absolutely  asymptotically  as 


(— )  for  large  Ji  . 
2.   Pg  >  0^ 

The  terms  of  series  (5-9)  behave  absolutely  asymptotically  as  — ttTo 
(—)   for  large  t . 

5-   ^B  =  ^C 

The   terms  of  series  (5-9)  behave  absolutely  asymptotically  as   ^ /. 

for  large  /£  . 
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Appendix  (l) 

A  general  Slater  function!-  J  denoted  by  \|r   ,    (r  5.^6  „,0  ~:k)  is 


given  by 


-kr 


C(n,  je,m;k)r°;^  e  ^^   p|^"l(co8  9,.^) 


m8 


H5'  ^^'"m5 


J 


where  k  is  a  positive  real  nrmiber;   n,  /£,  m  are  integers  satisfying 
0  <Ji<  n  and  0  <   |m|  <  Xi   ^  -  A,B,   or  C;   and  5  =  1  or  2. 


C(n,^,m;k)  =  (2k) 


-i- 


(2Ul)U-\m\): 
2rt(;e+|m|):(2n): 


ni/2 


and  ^(0)  =  cos  in0  for  m  >  0,  ^(0)  =  —  ,  and  ^(0)  =  -sin  m0  for  m  <  0. 
Ihrovighout  this  paper  we  restrict  ourselves  to  Slater  fvmctions 


with  n  =  1  or  2  only.  These  are  the  following  five  functions: 

♦.  ^  ..  (r  .,e  .,0  _;k)  =  (kV«)^'^^  e  ""^^ 

i>^   -  -.  (r  .,9  _,0  _;k)  =  (kV5n)^/^  r  .  e   ^^ 
2,0,0;u-    (io  no"^n6'     v  /^  /    ^g 


5/  ,1/2^   ^"%6 


*o  1  r^.  (r  R>©  R/0  R^ls:)  =  (kVn)^  r  _  e   "-^^  cos  9 


H5 


.5/^x1/2^ 


-kr 


U5 


to  T  I   (r  ;,,9  -,,0  _;k)  =  (kV")  '  r  ^  e   ^"  cos  9  -  cos  0  _ 


5/,^l/2. 


-kr 


H5„ 


to  T  T.  (r  ^,9  -,0  _;k)  =  {k^/nr'r  ^   e   ^  cos  9  -  sin  0  _  . 


These  are  denoted  by  t   (5),  for  convenience,  and  the  constants 


C(n,  x,m;k)  are  denoted  by^C. . 


The  expansion  for  r^  is  given  as  fellows'-  -' 
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/  -  ,  ^  /   £+1 


'll     =  £  r^  D^J-A<  /  ^A>  )  ^^^°^  ^Al)^(<^°^  ^A2)^°^  -(^A1-'^A2) 


where  D^^  =  1,  D^^  =  2  [j^   ,  m  ^  0;  r^<  =  min(r^,r^),  r^>  = 


"^^'■a1''*A2^- 

dv.,  5  =  1  or  2  is  given  as  follows: 

•^^5  =  ^A5  ^^°  ®A5  "^^AS  ^®A5  ^^Ab' 

[2l 
Ibe  P7(cos  e)  are  the  standard  Legendre  functions'-  ■• 
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Appendix  (2) 
The   fomrulas  used  to  expand  functions  with  arguments  (r„, ©„, 0_)  and 
{r„,Q„,0„)   in  terms  of  functions  with  arguments  (^a>®a>0a)  ^^®  "^^  following: 

03-  0, 
rg  cos  Gg  =  pg  -  r^  cos  9^ 


r-'  -  E  (p^  /  pr>)  ^i(-«  «a) 


n-1  -^B  ^'-°° 

r^   e 

®  -  "     i=0 


(Pb-a)''/'£  Eio^(^°^  ^a)  ^ni^^'V^B) 


r^  sin  Gp  sin  0^  =  r^  sin  9^  sin  0^ 


r_  sin  9_  cos  0_  =  r  (cos  a  sin  9.  cos  0.  +  sin  a  cos  9.) 

r_  cos  9_  =  Q„  -   r  (-sin  a  sin  9.  cos  0.  +  cos  a  cos  9.) 

^c^  =  EC  Eij(pJ<  /pJ>^)Pi(cos  a)P^(cos  9^)cos  J0 


J' 
^Va^'^^^  £  t  E^jpJCcosa)pJ(cos9^)co8  J0A^j,i(k, ^Pg) 


n-l  -^C         ''-«'  i 
^  -  "      1=0  J=0 


Here:     p     ^    «     mln(p  ,r,):  p  =     max(p  ,r.),  ji  =  B  or 


^io  °  ^'      Hi  '  2(i-j).'   /  (i+j);  J  ^  0 

E^^  =  2i+l;   E^j  -  2(2i+l)(i-j).'   /  (i+j).'     i  ^  0 


54 


^^i(k,r^,p)  = 


Sn 


a(-k)' 


(7j^(k,r^>p)) 


7^{^,r^,p)     =     I   i(kr^)K   ^^(kp)       ^a  -  ^ 


i+ 


i^2 


=   I    (kp)  K    (kr  )       r  >  p. 
i+  I      i+  2 

I   ^  and  K   .  are  the  standard  Bessel  fiinctions  of  Imaginary  argument 
and  half  integral  orders.  For  convenience,  ^^ .(k,r  ,p)  and  ^„.(k,r  ,p) 


li' 


'2i' 


are  denoted  by  p. (k,r.,p)  and  q.(k,r.,p)  respectively,  sometimes  written 
p.  or  q.  for  convenience. 

The  expansions  relating  functions  of  {t„,9„,0„)   to  those  of 
(r.,9.,0.)  are  easily  derived  from  those  relating  functions  of 
(r„,9^,0„)  to  those  of  (r  ,0.,0.).  Ihe  latter  expansions ^ '  J  easily 
follow  from  the  elementary  properties  of  Bessel  functions  and  Legendre 
polynomials . 

Some  relations^-!  between  Pj,q.  ,7^  and  ^t.  used  in  this  paper  are 
given  here. 


kr^p 


p^(k,r^,p)  ^^j^ 


7^.^(k,r^,p)  -  7i^i(k,r^,p) 


kr^p 
q^(k,r^,p)  =  2j^ 


Pi.]L^k,r^,p)  -  Pi+i(k,r^,p) 


-1 


-k"  p^(k,r^,p) 


2r^p 
C5i(k,r^,p)  -  21^ 


Pi.l(k,r^,p)-p^_i_^(k,r^,p) 


kr^ 
2i+l 


q^_^(k,r^,p)-q^^^(k,r^,p) 
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AppencLlx  (3) 

The  following  formulas  are  used  to  integrate  over  the  ©^^,  9^ 

variables  in  the  Integrals.  Itqrconsist  essentially  of  expansions  for 

sin  9  cos  9P^(cos  9)  N,  M,  n  non-negative  integers  in  terms  of  linear 
n 

ccanbinations  of  Legendre  polynomials,  and  fonmilas  for  expanding 

sin  9  P^{cos  9)  and  cos  9  P^(cos  9)  in  terms  of  Legendre  fvmctions 

N' 
P  (cos  9)  of  lower  upper  index  N'.  They  can  easily  be  derived  from 
n 

the  recurrence  relations  for  Legendre  functions  and  their  differential 
equations  ■-  -^ . 

The   general  expansion  formula  is 

sin  9  cos  9  P;(cos  9)  =  ^I  G^i  Pn-N-h*f2i^*=°^  ^^ 

where  the  G  ._..  (written  G  for  short)  are  rational  functions  of  n. 
The  first  13  G's  are  listed  here: 

G,_  »  ^^  G        »+l 


nlOO    2n+l  nlOl     2n+l 

_  n(n-H)  _  _  P(n-H) 

nOlO  ~   2n+l  nOll  "    2iH-l 

n(n-l)  2n^+2n-l  (m-l)(n+2) 

n200  "(2n-l)(2iH-l)        n201  ~  (2n-l)(2n+5)      n202  ~  (2n+l)(2n+3) 

r      n(n-l)(n+l)  p(ptl)       n  -n(n+l)(n+2) 

nllO  ~  (2n-l)(2n+l)       nlll  ~  (2n-l)(2n+3)     nll2   (2n+l)(2n+5) 

(n-l)n(n|l)(n-.-2)  „  (n-l)n(nH-l)(n-|-2)       ,  (n-l)n(n-fl)(n-f2) 
n020  ~    (2n-l)(2n+l)     n021  ~   (2n-l)(2iH-3)     n022  "  ( 2iH-l ) ( 2n+3 ) 

In  general: 

(Polynomial  in  n  of  degree  <  Mf2N) 

nMNi  ~  (Polynomial  in  n  of  degree  =  M+N) 
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In  particiilar: 

„  (Polynomial  In  n  of  degree  =  W-2N)   ,  _  q  ^j.  u+„ 

nMNj  "^  (Polynomial  in  n  of  degree  =  MfN) 

In  addition  we  shall  use  the  following  formulas: 

* 

sin  e  P***'''(cos  e)  +  2N  cos  9  Ocos  9) 
n  n 

+  (n-H+l)(n+N)sin  6  P^"-'-(cos  9)  =  0 

cos  e  I^(cos  9)  =  G^^  P^.^(cos  9)  +  G^q^I^^^(cos  9) 

-  N  sin  9  P^''''"(cos  9). 
n  ^ 


* 
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Appendix  {k) 

We  shall  derive  some  simple  Inequalities  for  the  Bessel  functions 

I   ^(x),  K   j^(x)na  positive  integer  and  x  any  resil  number  greater 

than  0.        Starting  with  the  integral  representation  of  I   ^(x)'-  -I 

1  """^  2 

I  ,{x)     =     —^ I  e"^^l-t2)°dt, 

n+ 


clearly 


i         v;?2^^/2n:   4 


n.1  1 

X  e  f  /,  ^2xn, 


yjt  2    /  n:  /^  n+  2 


^2  X    ,1 


7^2^ 


n: 


-1 

To  evaluate  J  (l-t^)°dt  we  write  the  integral  as  j  (l+t)°(l-t)°dt  and, 
-1  ^1 

by  repeated  integration  by  parts,  we  obtain 

1 


-1 


1 


Hence 


(^)         e  ^nl  ^  J    (^j  ^  (2x)       Vn! 
y«  (2n+l):       n+  I       yrt  (2n+l): 

Taking  the  integral  representation  of  K   ^(x)'-  •• 

^   l^''^  =     n-h  1/9  J  ^   ^*  -^^  ^* 
"**  2       n:  2       1 


58 


we  set  u+1     =  t  and  obtain 


(D  00 


j    e-^^t^-D'^dt     =     e-^l     e-^(u^+2u)du 


0 

00 

X   (       -X     2n,  (2n).'  ^-x 

e       u     du    =    ',.  (,   e 


>e        .     .       2n+l 

0 


On  the  other  hand 

OD  00  00 

.-xt^2n^^  _  i 

2n+l 


'       -xt,^2  ^ vn,^  ^   f   -xt^2n,^  ^  (       -xt^2n,.    (2n)l 
e   (t  -1)  dt  <     e   t  dt  <     e   t  dt  =  p— fr  ; 


J 

1 

so  we  have 

-x 


1  0  ^ 


(2")'.)^       <     K  (x)     <     IMliV- 

,n+  1/2    ,     -    ^  r^     -     /^xn+  1/2    , 
)  '    nl  n+  ;r  (2x)  '    n.' 


Taking  the  definition  of  y   (k,x,p),  p  >  0,  k  >  0  (Appendix  (2)) 
we  can  easily  derive: 


-k(x+p)     „  n+  -  kx  n+  2 

%^^kl)  <     7,(k,x,p)     <     l^(^)  forO<x<p 

-k(x-l-p)    n+  -  kp  ^  n+  - 

^-^r— r—  (^)   ^  <  7  (k,x,p)  <^(^)      for  x>p. 

2n+l   ^x  —   n^  '  *         —  2rH-l  ^x 


In  the  series  of  integrals  of  Sections  5A^5  and  6  we  replace  the 
7's  by  these  majorant  and  minorant  functions  to  obtain  majorant  and 
minorant  series  of  integrals.  ¥e  then  can  find  lower  and  upper  bounds 
for  the  terms  of  these  minorant  and  majorant  series,  respectively,  by 
the  following  inequalities: 


-  39  - 


-kx  n+1              /       ,                         n+1 
ex  .      I      -ky  n,         .    x 


ttfi 
"br  <- 


,  kx  n+1 

ky  n,       ^     e     X 


OD 


e'Vdy    < 


-kx 


(n-l)x 


n-1 


n  a  positive  Integer  greater  than  1,  and  x,k  real  rnmbers  greater  than  0. 
We  shall  also  use  the  following  definitions  in  this  paper: 


r+ 

^ 

x; 

X     > 

0 

b^ 

— 

0; 

X     < 

0 

x_ 

= 

x; 

X      < 

0 

X 

s 

0; 

X       > 

0 
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